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ABSTRACT 



We find evidence, by both observation and analysis, that crater ejecta play 
an important role in the small crater (less than a few km) populations on the 
Saturnian satellites, and more broadly, on cratered surfaces throughout the 
Solar System. We measure crater populations in Cassini images of Enceladus, 
Rhea, and Mimas, focusing on image data with scales less than 500 m/pixel. 



We use recent updates to crater scaling laws and their constants (Housen and 



Holsapple, 2011[ ) to estimate the amount of mass ejected in three different 



velocity ranges: (i) greater than escape velocity, (ii) less than escape velocity 
and faster than the minimum velocity required to make a secondary crater 
("Wmin), and (iii), velocities less than fmin- Although the vast majority of mass 
on each satellite is ejected at speeds less than fmin, our calculations demon- 
strate that the differences in mass available in the other two categories should 
lead to observable differences in the small crater populations; the predictions 
are borne out by the measurements we have made to date. In particular, 
Rhea, Tethys, and Dione have sufficient surface gravities to retain ejecta mov- 
ing fast enough to make secondary crater populations. The smaller satellites, 
such as Enceladus but especially Mimas, are expected to have little or no 
traditional secondary populations because their escape velocities are near the 
threshold velocity necessary to make a secondary crater. Our work clarifies 
why the Galilean satellites have extensive secondary crater populations rela- 
tive to the Saturnian satellites. The presence, extent, and sizes of sesquinary 
craters (craters formed by ejecta that escape into temporary orbits around 



Saturn before re-impacting the surface ( Dobrovolskis and Lissauer, 2004; Al- 



varellos et al. , 2005 Zahnle et al. , 2008 ) is not yet well understood. Finally, 



our work provides further evidence for a "shallow" size-frequency distribution 
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(slope index of ~ 2 for a differential power-law) for comets a few kilometers 
diameter and smaller. 



Keywords: crater ejecta, secondary craters, sesquinary craters, Saturnian 
satellites 
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1 Introduction 



Without access to planetary samples from outer Solar System objects, cata- 
loguing the size-frequency distribution (SFD) and spatial distribution of im- 
pact craters is a common means to determine surface ages. Constraining age 
on an object, and between objects, is critical to establish a coherent and con- 
sistent chronology of events within a planetary system, or across the Solar 
System. Combining the measured crater SFD with known or modeled crater 
formation rates provides a means to estimate absolute age. 



To derive accurate ages using impact craters, in either the relative or absolute 
sense, one must first determine the sources of impactors that make craters. 
Impact craters can be primary, secondary or sesquinary; see Figure 1. Primary 
craters are made by direct impact of comets or asteroids. Secondary craters 
are the result of essentially ballistic trajectories of ejecta from the primary 
crater to some distance away. For typical impact speeds of hehocentric comets 



onto Saturnian satellites of several to 30 km/s (Zahnle et al. , 2003 Dones 



et al. , 2009), ejecta can be launched at speeds from a few hundred m/s to 
several km/s. When an ejectum is launched at a speed faster than the escape 
velocity of a moon, it can go into orbit about the planet. Most escaped ejecta 
are eventually swept up by the source moon, but the orbits of some escaped 
ejecta can be sufficiently perturbed, or the original ejection velocity so high. 



that the ejecta will impact another satellite (Alvarellos et al. , 2005, 2002). In 



either case, craters formed by ejecta that initially escape their parent object 



are called sesquinary ("l2-ary"; formerly "poltorary") craters ( Dobrovolskis 



and Lissauer, 20041 Zahnle et al. , 2008). Because secondary and sesquinary 



craters are products of primary craters, and because the larger (and therefore 
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generally older) primary craters create the most ejecta, older terrains will have 
the greatest number of craters of all types. This introduces uncertainty in the 
number of primary craters, which is the only kind to be trusted chronometri- 



cally (McEwen and Bierhaus, 2006) 



Most Voyager-era studies of craters on the Saturnian moons identified two 



basic classes, called Population I and Population II (Smith et al. , 1981, 1982) 



Population II follows a steeper size-frequency distribution than Population I 
(i.e., in relative terms. Population II has a smaller number of large craters). 



One hypothesis proposed (e.g.. Chapman and McKinnon (1986)) that helio- 
centric comets are responsible for Population I, while planetocentric debris 
from cratering or catastrophic disruption of moons is responsible for Popu- 



lation II. A second hypothesis (e.g., Strom and Woronow (1982)) proposed 



that planetocentric debris largely is responsible for both populations. A third 



(Plescia and Boyce, 1985) suggested that Population II was created by late- 



arriving bodies from heliocentric orbit. The existence of two populations was 
questioned by Lissauer et al. (1988), who argued that the Voyager data on 
small craters might be consistent with (near-)saturation equilibrium produced 
by heliocentric impactors (Hartmann [1984 ). For our purposes, the key point 
is that ever since craters were first seen on the Saturnian moons, the question 
of whether the impactors were heliocentric, planetocentric, or both has been 
controversial. 

Since the Voyager fiybys, a number of developments have taken place (see 



reviews by Chapman and McKinnon (1986); Schenk and Pappalardo (2004) 



Dones et al. (2009)). First is the confirmation of the Shoemaker and Wolfe 



(1981) and Shoemaker and Wolfe (1982) prediction that heliocentric "ecliptic" 



comets should dominate primary impacts on planets and regular satellites in 
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the outer Solar System (Levison and Duncan 1997 Zahnle et al. 1998 2003[ ), 
and the discovery of their putative source region, the Kuiper Belt/Scattered 



Disk (Barucci et ah, 2008). Second is the resurrection of Shoemaker's (Shoe 



maker, 1965) observation that secondary craters can be a major component 



of the small crater population (Bierhaus et al. 2005 McEwen et al. 2005 



McEwen and Bierhaus 2006). Third is the detailed modeling of the dynami- 



cal behavior of escaped ejecta in multi-body systems (such as the Saturnian 
system), leading to the first quantitative estimate of the effect of sesquinary 



ejecta launched from large craters on the mid-sized moons (Moore et al. , 2004) 



on the cratering population (Alvarellos et al. , 2005). Fourth is the ongoing lab 



oratory and analytical work into the physics of crater ejecta (Melosh, 1984 



Housen and Holsapple 2011). Fifth is the possible origin of irregular satellites 



of the giant planets, and possible consequences of their coUisional evolution 



for the cratering record on inner satellites (Bottke et al. , 2010). The cap 



stone, of course, is the presence of the Cassini mission at Saturn, and the 
imaging (ISS) data that enable a significant improvement of the catalogue of 
craters. Crater counts on the Saturnian moons using ISS images have been 



published for the following moons: Phoebe and lapetus (Porco et al. , 2005) 



Hyperion and Phoebe (Thomas et al. , 2007); basins on the leading face of la- 



petus (Giese et al. 2008); Enceladus (Porco et al. 2006 Kirchoff and Schenk 



2009); and Mimas, Tethys, Dione, Rhea, lapetus, and Phoebe (Kirchoff and 



Schenk, 2010). Counts for Titan, using Cassini's radar instrument, have been 



published by Lorenz et al. (2007); Wood et al. (2010), and Neish and Lorenz 



(2011). 



The first development cited above - the prediction that ecliptic comets are the 
primary impactors on the regular satellites of the giant planets - is based, in 
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large part, on the sheer numbers of comets that seem to wander the outer So- 



lar System (Zahnle et al. 1998, 2003). Estimates of the impactor populations 



at each planet mostly rely on interpolation between observations of Jupiter- 
family comets (generally seen well interior to Jupiter's orbit, and prone to 
disintegrate for reasons that are not understood) and counts of (large) Cen- 
taurs and Kuiper Belt Objects, corrected for discovery biases with the help of 



dynamical models (Dones et al. , 2009). Detailed comparison of crater SFDs 



and the independently known population of impactors, such as has been car- 



ried out for the Moon by Marchi et al. (2009), is not possible for satellites in 
the outer Solar System. 



On synchronously rotating satellites, ecliptic comets are much more likely to 



strike the moons' leading hemispheres (Shoemaker and Wolfe, 1982; Horedt 



and Neukum, 1984a; Zahnle et al. , 2001). The expected "apex-antapex asym- 



metries" do not occur to the extent expected, although craters on bright 



terrains on Ganymede (Zahnle et al. 2001 Schenk and Pappalardo 2004) 



and rayed craters on Saturnian satellites (Schenk and Murphy, 2011) do ex- 



hibit smaller asymmetries in the expected sense. (Triton is just weird; see 



Schenk and Zahnle (2007) if you must know the details.) The absent or muted 



apex/antapex effect might be due to crater saturation, true polar wander, or 
some form of nonsynchronous rotation. In addition, moons closer to their par- 
ent planets should be more heavily cratered because of gravitational focusing 



(Smith et al. , 1981 , 1982; Zahnle et al. , 2003). This effect is not observed either 



perhaps as a result of saturation, at least for small craters, due to an early era 



of heavy bombardment (Chapman and McKinnon 1986 Dones et al. 2009 



Richardson, 2009) 



Alternatively, one can abandon the hypothesis that the impactors have a he- 
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liocentric origin, and posit that the craters are produced by planetocentric 



(Saturn-orbiting) bodies (Horedt and Neukum, 1984b Neukum, 1985). The 
planetocentric model has the advantage that planetocentric cratering only 
weakly favors a satellite's leading or trailing side, depending on whether the 



debris fall from outside or inside the moon's orbit, respectively (Horedt and 



Neukum, 1984a Alvarellos et al. , 2005). However, planetocentric bodies on 



crossing orbits are generally short-lived (Burns and Gladman, 1998; Alvarellos 



et al. , 2005), so they must be resupplied if they are to contribute to latter-day 



cratering. Ironically, the most plausible source of resupply is impact by helio- 
centric comets, which might produce copious fragments that rain back onto 
the moons. This realization motivates our study of the role of ejecta from 
cometary impacts in producing the observed crater populations. 

The outline of our paper is as follows. In Section 2, we apply the cratering 



ejecta model of Housen and Holsapple (2011 ) to derive the expected mass avail- 



able to form sesquinary and secondary craters on Saturn's mid-sized moons 
and Jupiter's three icy Galilean satellites. In Section 3, we describe the image 
processing we performed on the publicly available Cassini images we used to 
measure crater populations. In Section 4, we present the size-frequency distri- 
butions of craters on terrains on Enceladus, Mimas, and Rhea. In Section 5, 
we compare the observed and predicted populations of sesquinary craters on 
the Saturnian moons. Section 6 summarizes our conclusions. 



2 Estimating the Role of Ejecta as Impactors 

The crater formation and evolution process has been long approximated by 
a set of scaling laws, anchored to physical reality by laboratory-scale exper- 



9 



iments, nuclear explosions, observations of impact craters across the Solar 
System, and increasing sophistication of numerical modeling. Although the 
details of any one impact are not captured by the scaling laws, they provide 
a reasonable description of trends. 



For the current research, we focus on the mass-velocity relationship of crater 
ejecta. In particular, observations of laboratory-scale cratering experiments (?) 
reveal that there is an inverse relationship between ejecta mass and ejection 
velocity, which is to say that less mass is ejected at higher velocities, see 
Figure 2. Early stage ejecta from near the point of impact move the fastest, 
but are the smallest portion of ejected mass. Late-stage ejecta are moving the 
slowest, but are the bulk of ejecta mass. Drilling at terrestrial craters and 
observations of craters on other surfaces indicate that the end of the crater 
excavation phase is a transition to material that isn't so much ejected as it is 
an overturned flap of the surface. 



Using the most recent summary of laboratory measurements of cratering ejecta 



from Housen and Holsapple (2011), and impact velocities for ecliptic comets 



on the satellites from Zahnle et al. (2003), we estimate the total and fractional 
ejecta masses for a 1 km diameter cometary (density of 600 kg/m^) impactor 
on various Saturnian satellites. We use a 1 km comet because it is sufficiently 
large to generate ejecta masses that contribute to the observable secondary 
and sesquinary crater populations, but small enough that such impacts should 
occur on geologically short timescales. Table 1 lists properties of the satellites 
used in our calculations. To complete the calculations, we must quantitatively 
determine the appropriate velocity regimes that separate the ejecta blanket, 
secondaries, and sesquinaries. We discuss these velocities next. 
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2.1 Making an Impact: v, 



We define Vmm as the boundary between the ejecta blanket and secondaries, 
i.e., it is the minimum velocity at which an ejectum can make a secondary 
crater. 

We use observations of adjacent secondaries on Europa, and the simple planar 
baUistics equation, to set plausible values of for icy surfaces. The baUistics 
equation is: 



where r is range, v is ejection velocity, 6 is the ejection angle (measured relative 
to local horizontal), and g is the surface gravity. Measuring r enables us to 
calculate v: 



First, we examine Rhiannon, a 15 km diameter crater on Europa; see Figure 3. 
Rhiannon is the smallest primary on Europa imaged with sufficient resolution 
to see its adjacent secondary crater field. (Galileo imaged smaller primaries on 
Europa, but not with sufficient resolution to determine whether or not they 
have secondaries.) The point from which the fragments that make these secon- 
daries originate is somewhat uncertain. The late ejecta come from the region 
close to the crater rim. The final crater diameter is larger than the transient 
diameter due to collapse of the crater after the excavation phase is complete. 
We use a point roughly halfway between the crater center and the final crater 
rim as a reasonable location for the origin of the secondary fragments. The 



sin 29 
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(2) 
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distance between our presumed origin point and the first appearance of tlie 
adjacent secondary crater population is about 17 km. Substituting r = 17 km, 
9 = 45°, and g = 1.31 m/s^ for Europa into Equation 2 gives v ~ 150 m/s. 

Next, we examine Tyre, tfie largest known primary impact structure on Eu- 
ropa; see Figure 4. Depending on the location of the crater rim. Tyre is 
45 — 50 km diameter. Following the procedure for measuring range we defined 
for Rhiannon, we find that secondaries begin to appear about 50 km from the 
estimated point of origin. Using 50 km in Equation 2 gives v ~ 250 m/s. 

Thus, in our subsequent analysis, we examine two different cases for the min- 
imum velocity required to make secondaries, Vmm = 150 m/s and fmin = 
250 m/s. Although the actual value could differ for the Saturnian satellites, 
the general trends described here would not change. 



2.2 The Great Escape: v^^ 



Alvarellos et al. (2002, 2005) clarified that the relevant velocity to escape a 



moon in orbit around a planet isn't the classical escape velocity appropriate 
for an isolated body, but rather the velocity required to reach the moon's Hill 
Radius, which is: 

1/3 



Rh — O-r 



Mrr 



3{Mm + Mp) 



(3) 



where is the moon's semi-major axis. Mm is the mass of the moon, and 
Mp is the mass of the planet. The velocity necessary to reach Rh is then: 



Rjj — RuRm 

. Rl - R'J^ 



1/2 



(4) 
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where Rm is the moon's radius, fesc = \J'^GMm/ Rm is the classical escape 
velocity, and is the moon's mass. When estimating the amount of mass 
that escapes a moon, we use rather than fesc- In general, v^^^ < v^sc- 



2.3 Scaling Law Calculations 



Our calculations are based on the expressions and definitions within [Housen 



and Holsapple (2011). First is the expression for the cumulative mass ejected 



faster than a given velocity for a crater in the gravity regime, defined in terms 
of impactor properties: 



M{v > v) 
rrii 



Vi) \Pi) 



(5) 



where C4, z/, and n are constants, and M{v > v) is the mass ejected faster 
than velocity v, Vi is the impactor velocity, pt is the density of the target's 



crust, and pi is the impactor density. Housen and Holsapple (2011) indicate 



that u is about 0.4 regardless of target type. The value of the exponent p 
is expected to lie between ^("momentum scaling") and | ("energy scaling"). 



For the materials tabulated by Housen and Holsapple (2011), p ranges be 



tween 0.35 for perlite/sand mixture to 0.55 for water and rock. The values 
of p and C4 for ice are not given, but K. Housen (personal communication, 
2011) recommends using the values for weakly cemented basalt as the closest 
approximation for ice {p = 0.46, C4 = 6.72 x 10~'^). We set pt = 900 kg/m^, 
appropriate for nonporous ice. 

For a 1 km diameter comet with a density of 600 kg/m^, the impactor's mass 
is 3.14 X 10^^ kg. We use Zahnle et al. (2003) as the source for typical values 



of Vi for ecliptic comets striking the different satellites. Substituting numer- 
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ical values into Equation 5, we obtain the following for ejecta masses on icy 
satellites: 



/ ^ \ -1.38 / \ -0.2 

M{v >v)= 6.72 X IQ-^rm l^j Ijj (6) 



We will henceforth write M{v > v) as M{v). First we calculate the ejected 
mass that escapes the moon, M{v^^^). In our discussion, we assume that this 
mass is equivalent to the amount of mass available to make sesquinary craters, 

Mi.5. 

/ H \ -1-38 / s -0.2 

Mi.5 = M{v^J = 6.72 X lO-'m, f ^ j i^ j (7) 



Alvarellos et al. (2005) demonstrated that, for the mid-sized Saturnian satel- 
lites, most but not all of the sesquinary fragments re- impact their source moon. 
Their results were bracketed by 81.5% re-impact (ejecta from Odysseus crater 
on Tethys) on the low end, and by 99.6% re-impact (ejecta from Herschel 
crater on Mimas) on the high end. 

Next we calculate the mass available to make secondary craters. Msec, which 
is: 

Msec = M{v^^) - Mi,5 (8) 



where M{v-[ain) is the mass ejected faster than the minimum speed necessary to 
make secondary craters, fmin- This value is likely variable, dependent on target 
properties, but observations of Europan secondaries provide some constraints 
on the minimum possible value. See our discussion in Section 2.1, where we 
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identify two plausible values for Vmin of 150 m/s and 250 m/s. 

— 1 38 / \ 

M(^^,„) = 6.72 X 10-=^m, (^) ' 1^ ] (9) 



Later we find it useful to compare these masses with the total ejecta mass 



for a crater, Mtot- As Housen and Holsapple (2011) discuss, the volume of the 



final transient crater - at the end of the crater excavation phase but before 
modification due to wall collapse, etc. - has a corresponding mass: 

-^crater = ^cratcrPi-R^ (10) 



where kcrater is a constant, and R is the radius of the transient crater at 
the end of the excavation stage. This value does not correspond to an ac- 
tual physical mass, since some significant fraction of any crater is formed by 
compaction and/or inelastic compression rather than excavation. K. Housen 
(personal communication, 2011) recommends that Mtot/Mcrater is about 0.5 
for craters on icy satellites. Thus, we use the following expression for the total 
mass ejected: 

Mtot = 0.5 Merater = 0.5 /Cerater Pt R\ (H) 



where we assume krratpr = 0.6. 



To calculate R, we used the expression from Housen and Holsapple (2011) 
for the diameter of the transient crater in the gravity regime, given a certain 
impactor: 

2-\-fi — Gu _ jj, 

D = 2Hj£Ly"(ei] '''"" fs|V"', (12) 
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where Hi is a constant. K. Housen (personal communication, 2011) recom- 
mends using Housen and Holsapple (2011 )'s value for sand of Hi = 0.59 



because of the similarity in porosity, although there may be differences due to 
friction angle. Substituting fi = 0.46, u = 0.4, we have 



/ 2\0-iS7 / s 0.325 
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(13) 



as compared with 



D = 1.1 



' 2\ 0.217 / .\ 0.333 

f J \ I Pi COS t \ 



0.783 



Pt 



(14) 



the expression used by Zahnle et al. (2003); Alvarellos et al. (2008); Zahnle 



et al. (2008), where i is the angle of incidence (0 for a vertical impact). These 
authors assume /i = 0.55, appropriate for impacts into wet sand, rock, and 



water (Schmidt and Housen 1987 Housen and Holsapple 2011). We incor- 



porate impact angle into our calculations (Equations 6 and 13) by assuming 
an average primary impact angle of 45°, and scaling the impact velocities in 



Zahnle et al. (2003) by cos 45° = 0.707. 



With Mtot, it is a simple exercise to calculate the mass available to make an 
ejecta blanket (Mbik): 



Mbik = Mtot - M{v^i^) 



(15) 



Figure 2 is a schematic illustration of these relationships. Any mass moving 
faster than t>^^ is available to make sesquinary craters. Mass ejected faster 
than fmin but slower than makes secondary craters {Msec)- Mass ejected 
slower than t^mm is available to make the ejecta blanket (Mbik). The smaller the 
difference between Wmin and w^^,, the less mass is available to make secondary 
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craters. (For very small moons, we have fmin < v^^, in which case no mass is 
available to make secondaries.) Conversely, the greater the difference between 
the two values, the more mass is available to make secondary craters. 

The effective escape velocity v^^^ is a function of the mass and radius (or bulk 
density and radius) of the target object, as well as the moon's distance from 
its planet (see Eqs. 3 and 4). The minimum velocity of ejecta that can make 
secondaries, fmin, is a function of the surface strength. Presumably is 
greater for rocky objects than for icy objects, although the difference may not 
be significant. For the mid-sized icy satellites of Saturn, fesc is smaller than it 
is for the terrestrial planets (or the Galilean satellites). 

For purposes of comparison, we calculate Mtot, M^k, Msec, and Mi. 5 for a 1 km 
diameter comet impact on several mid-sized icy Saturnian satellites and the 



three icy Galilean satellites. We use the impact velocities from Zahnle et al. 
(2003), and an impactor density of 600 kg/m^. 



In all cases, the vast majority of the ejected mass resides in the "ejecta blan- 
ket" . Table 2 lists the absolute mass and fractional mass for each moon when 
Vmin = 150 m/s, the case which minimizes the amount of mass available for 
Mbik. Even in this case, between 95% and about 99% of the mass is in Mbik- 
However, enough ejecta are launched at sufficient speeds that such an impact 
should add to the observable secondary and sesquinary crater populations. To 
illustrate this point further. Table 3 lists Msec and Mi, 5 in terms of the im- 
pactor mass. Depending on the value of f for each moon, the primary impact 
can eject multiple impactor masses to create secondaries and/or sesquinaries. 

Figures 5 and 6 plot the fractional and absolute mass available to make sec- 
ondary craters on each of the satellites. Mimas has no mass available to make 
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secondary craters in either case, because v^^ < Vram for that moon, while the 
same is true for Enceladus in the Vmm = 250 m/s case. In both cases, Europa 
has the highest mass available for secondary craters. Tethys, Dione, and Rhea 
are the most similar to one another among the satellites addressed in this 
study. 

Figure 7 plots the fractional and absolute mass available to make sesquinary 
craters on each of the satellites. Mi. 5 is independent of fmin, and so is the same 
between the two values of fmin used here. Mimas and Enceladus have orders 
of magnitude more mass available in this ejecta category than the Galilean 
satellites, and have significantly greater amounts of sesquinary ejecta than the 
other mid-sized Saturnian satellites. 



2.4 Secondaries: Near and Far 



Frequently secondaries are grouped into two basic populations: adjacent and 



distant (McEwen and Bierhaus, 2006). Adjacent secondaries are those that 
appear immediately outside the continuous ejecta blanket of a primary crater, 
and are a distinct high-spatial density annulus of small craters around the 
primary. The "distant" secondary is a broad classification, covering those sec- 
ondaries that appear several crater radii away, to those that are hundreds (or 
even thousands) of km away (McEwen et al. , 2005; McEwen and Bierhaus| 
2006] [Prebhch et al.\ [20071 |Dundas and McEwenl [20071 |Robbins and Hynekj 



2011). Of course, an ejecta fragment with a given velocity will result in dif- 



ferent ballistic ranges on objects with different surface gravities and sizes. 
Formally, the ballistic range R of an ejecta fragment, launched from a point 
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on a sphere and traveling less than the escape speed of the object, is: 



R = 2(j)Rm (16) 

where Rm is the radius of the moon, and is the half-angle distance of travel, 
defined by: 

vl^ sin 6 cos 6 
gRm - f cj cos^ 6 



(Vickery, 1986), where Wej is the fragment ejection velocity, 9 is the ejection 



angle, and g is the surface gravity of the moon. 

An illustrative comparison is to normalize the range by the circumference of 
the body, i.e. R/Cm, where Cm = '^T^Rm is the circumference of a moon. In 
Figure 8, we plot R and RjCm for t>cj = 150 m/s, i.e., the lower limit to fmin, 
and in Figure 9 we plot the same information for fgj = 250 m/s, the maximum 
value of Wmin we use for the current analysis. 

For fmin = 150 m/s, a secondary crater on Enceladus will form no closer 
than over 300 km distant, and even for the more massive Rhea, the closest 
secondary forms over 80 km away. (In this model, no traditional secondaries 
form on Mimas for v-rmn = 150 m/s.) This minimum range is a few percent of 
Rhea's circumference and almost 20% of Enceladus's circumference. 

For fjnin = 250 m/s, the closest secondaries move further away from their 
parent primary. (For v^^m = 250 m/s, the available mass to make secondaries 
is zero for both Mimas and Enceladus.) The closest secondaries on Rhea are 
over 200 km distant, while those on Dione and lapetus are over 300 km distant, 
and those on Tethys are over 600 km distant - 20% the circumference of the 
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moon. 

On objects with small surface gravities, then, there may not be "adjacent" 
secondary populations in the traditional sense. Certainly we don't expect to 
see dense, overlapping fields of secondary craters around a large primary. This 
is true for the Saturnian satellites under discussion here, as well as asteroids, 
Kuiper Belt Objects, and other bodies where Vmin is close to, but still less 
than, Vesc (or J. 

2. 5 Discussion 

Based on examination of the tables and figures, we divide the satellites into 
different groups based on the relative values of Msec and Mi. 5. (Table 3 divides 
the satellites into groups with horizontal lines to illustrate the discussion.) In 
particular: 

Mimas and Enceladus: These moons have the weakest surface gravities and the 
smallest escape velocities of those under consideration, as well as the highest 
cometary impact velocities. This combination means that an impactor will 
create the largest crater on these satellites, and further, will generate more 
mass for sesquinary craters than for secondaries. 

Tethys, Dione, and Rhea: Larger than Mimas and Enceladus, these moons 
have higher surface gravities and escape velocities. In addition, impact speeds 
are lower. On these objects, there's generally more mass to create secondary 

craters than sesquinary craters. 
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lapetus: Despite its similarity in size to Dione and Rhea, lapetus does not 
have much ejecta available for either secondaries or sesquinaries, due to the 
low primary impact velocity of around 6 km/s. The low impact velocity has 
two effects: the resulting primary crater and its ejecta mass are smaller, and 
a slower-moving primary impact generates slower moving ejecta. 

The Galilean Satellites: These objects are clearly distinct from the Saturnian 
satellites, in that they have higher fractions of mass available for secondary 
craters, and lower fractions of mass available for sesquinary craters. Their 
higher escape velocities mean they retain more ejecta as secondary craters. 
For Europa, the combination of high impact velocity (which generates a large 
primary with lots of fast moving ejecta) and higher surface gravity means 
there is a significant amount of mass available to make secondary craters. 
This explains why clusters of secondary craters are everywhere on Europa 



(Bierhaus et al. , 2005) - even though that moon has a relatively low density 
of primary craters - and yet clusters are rare to absent on the much more 
heavily cratered Saturnian satellites. 



3 Cassini Image Data 

We now turn to our analysis of Cassini images of three of Saturn's mid-sized 
satellites and our techniques for measuring craters on the moons. 

3. 1 Description of Image Processing 

We obtained the image data from the NASA Planetary Data System (PDS) 
website (http : // img . pds . nasa . gov/ ) , and performed all image processing 
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using the USGS planetary image software package ISIS version 3, which pro- 
vides several routines to manage Cassini ISS data. We used a combination 
of the Cassini-ISS-specific routines and additional ISIS routines to generate 
image mosaics that were the basis for our measurements. The processing can 
be considered to consist of two phases: the first phase is "image preparation" , 
which is a sequence of steps that generates a properly calibrated image; the sec- 
ond phase is "mosaic preparation" , which collates information from each cal- 
ibrated image to generate the mosaic. [For more information on each routine, 



see the ISIS website (http: //isis . astrogeology .usgs .gov/Application/ ) 



The following is a brief summary of the processing steps for the image prepa- 
ration phase. 

(1) ciss2isis: imports a Cassini ISS image into ISIS by collating the separate 
data file (.IMG) and label file (.LBL) from PDS into a single ISIS-format 
cube file (.cub). This is a Cassini-specific routine. 

(2) cisscal: performs radiometric corrections to the images. This is a Cassini- 



specific routine (West et al. , 2010) 



(3) spiceinit: attaches geometry information for the spacecraft, instrument, 
and target body to determine the geometric properties (e.g., resolution, 
latitude, longitude) of each pixel in the image. 

(4) trim: removes a user-specified number of pixels from a user-specified edge 
of the image. Because ISS images often have a one or two pixel border of 
invalid data around an image, we removed a two-pixel wide border from 
each edge of an image. 

(5) lowpass: applies a N x M moving filter to the image. ISS image data 
are often returned using a compression scheme that assigns a fixed data 
volume per two-line pair. The data volume is more than sufficient for a 
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single line, but not enough for the entire second line. The completeness 
of the second hne depends on the scene content (and thus the ability of 
the algorithm to compress the data), but is never total. This results in 
horizontal variably sized "jail bars", i.e., black bars with no image data, 
every other line in the image. To eliminate these data gaps, wc used a 
1x3 (line by sample) moving filter to generate image data. In other words, 
each data-gap pixel became an image pixel, with an intensity value that 
is the average of the pixels above and below it. See Figure 10. 
(6) camstats: outputs a text file containing summary information on the lat- 
itude, longitude, resolution, phase angle, and several other image prop- 
erties. While camstats does no image processing, it is very useful to have 
a text-searchable summary for each image. 

Once the previous steps have been applied to a set of images that comprise a 
mosaic, one must collect information on each image to determine the properties 
of the mosaic. The following is a summary of the routines and other steps used 
for the mosaic preparation phase. 

(1) mosrange: computes and outputs the latitude and longitude extent, and 
pixel resolution, of a set of images. The results are output to a ".map" 
file. The .map file contains the required latitude and longitude extent of 
the mosaic, the mosaic projection type, and other details necessary to 
reproject and combine the individual images into the mosaic. 

(2) qnet creates and edits a control network. This step requires a significant 
amount of user interaction. Briefiy, the user loads the set of images for 
the mosaic into qnet, which is a graphical user interface that allows the 
user to identify "tie points" between images. A tie point is a common fea- 
ture between two (or more) images. Although the geometry information 
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loaded into each image during the spiceinit step (mentioned above) con- 
tains roughly correct geometry information to align neighboring images, 
there are still residual errors at the pixel-level scale. The qnet routine en- 
ables one to identify the residual misalignment between the SPICE files 
(accessed using spiceinit) and the actual image alignment. An output of 
the user's interaction with qnet is a network file (a ".net" file), which 
contains the corrections to the image alignments. 

(3) jigsaw: more precisely aligns the images to sub-pixel levels. The inputs 
to jigsaw include the set of images to be mosaiced, and the network file 
created using qnet. 

(4) cam2map: using the data from the ".map" file, this routine converts each 
image from the 2D plane of the image detector to a map-projected image. 
ISIS provides several map projection options; we typically used sinusoidal 
(Mercator equal-area) projections to generate mosaics for crater measure- 
ment. 

(5) automos: creates a mosaic from a set of map-projected images. The output 
of automos is the product we used for crater measurement. 

Some image mosaics consist of images acquired from rapidly varying viewing 
geometry (i.e., during a close flyby), or during an extended viewing sequence. 
For these mosaics, the native image resolution and/or phase angle (important 
for identifying craters, or distinguishing them from other features) can vary by 
many tens of percent, or even more than a factor of two, between the first and 
last images of the mosaic sequence. In these cases, we do not use all images 
that belong to the complete sequence; instead, we use those images whose 
resolutions differ by no more than 50%. We make this choice because the ISIS 
routine automos creates an image mosaic with a single resolution (specified 
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by the .map file generated by the mosrange routine); it does not preserve the 
native resolution of the individual images. Minimizing the resolution differ- 
ence between images used to generate a mosaic ensures a relatively common 
completeness limit across the mosaic. 



4 Measured Crater Size-Frequency Distributions 



We display the crater size-frequency distributions (SFDs) as "R-plots" [rela- 



tive plots]. (Crater Analysis Techniques Working Group et al. , 1979) Recall 



that SFDs can often be approximated as power laws, or a series of power laws, 
such that dN/dD = kD^"^, where dN is the number of craters per unit area 
with diameters between D and D + dD, and k and q are positive constants. 
For example, primary craters on the terrestrial planets follow distributions 
with g ~ 3. To enhance structure, the R-plot divides dN/dD by a power law 
with g = 3, i.e.: 



_dN^ 
" ~ dD D-3 



At diameters less than a few km, terrestrial planet SFDs are roughly horizontal 
lines in R-plots, which is to say, g ~ 3. For subsequent discussion, we use the 
following terminology: if R increases with increasing diameter, the SFD is 
shallow, i.e., dominated by large craters. For example, if i? oc -D, dN/dD oc 
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4-1 Enceladus 



We focused our initial measurements on the young terrain of Enceladus, be- 
cause those regions should be the best representation of the true primary crater 
population, and the least affected by secondary and/or sesquinary craters. We 
also measured some regions of heavily cratered terrain to compare the crater 
populations between young and old(er) regions of Enceladus. Table 4 summa- 
rizes the regions we've measured, and discussion of the measurements follows. 

Young Terrain: We measured the crater populations within several regions of 
the young terrains of Enceladus, including portions of the 
ISS_003EN_LIMTOP004_PRIME, ISS_004EN_REGEO002_PRIME, 
ISS 01 1EN MORPH002 PRIME, and ISS 01 1EN N9COL001 PRIME mosaics 
(see Table 4). These mosaics span significant areas at and near the south pole, 
and a variety of longitudes (see Figure 11). 

The SFDs on young terrains on Enceladus show a few notable features (see 
Figures 13, 15, 17, 19, and 21). First, they all display a differential slope of 
roughly -2 (g ~ 2), regardless of location. Second, the crater density roughly 
correlates with latitude: the southernmost regions have the lowest crater densi- 
ties, while the lower latitudes have higher crater densities. Third is the "noise" 
in the crater SFDs; although the average trend for each region is a -2 differ- 
ential slope, there are density excursions ("bumps and wiggles") that depart 
from the linear trend. We have not yet ascertained what causes these depar- 
tures, but an initial explanation is that these density variations are evidence 
for a dispersed and generally unclustered secondary crater population. 

Old(er) Terrain: There are two key features of the crater SFD measured on 
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the older, more heavily cratered terrain. First, at diameters larger than a few 
km, the SFD has a -3 differential slope. Second, at diameters less than a few 
km, there is a clear transition in the SFD to a roughly -2 differential slope 
that parallels the SFD seen on the young terrains of Enceladus. 

The absence of a steeply-sloped SFD (i.e., g > 3) at small crater diameters, 
even within heavily cratered regions, supports our prediction that secondaries 
should have minimal influence on the small crater populations of Enceladus. 
The presence of a steeper SFD at larger diameters is an intriguing evolution 
of what we preliminarily predict to be a -2 differential production population. 
We raise two points to address this evolution. 

First, if any secondaries form on Enceladus, we expect they will be only the 
slowest moving, and therefore the largest. (Just as there exists a clear inverse 
mass-velocity relationship for the crater ejection process as a whole, there 
is also a more granular inverse mass-velocity relationship for the fragments 
themselves: higher-velocity fragments are smaller, while slower-velocity frag- 
ments are larger. This is why secondary craters trend to smaller diameters 
with increasing distance from their primary crater, even though the fragment 



ejection/impact velocity is growing. (See McEwen and Bierhaus (2006) for 



more discussion on this in regards to secondaries, or Bart and Melosh (2010) 
in regards to ejected boulders.) In addition, these slow moving fragments will 
not re-impact near their parent primary, as discussed in Section 2.4. Rather, 
they will travel some significant distance away, and may be unrecognizable 
as secondaries due to lack of the distinct morphologies of tightly packed sec- 
ondary craters. The end result is a population of craters that: (i) follow a 
steeper SFD but are not clearly secondaries surrounding a certain primary; 
and (ii) extend down to some minimum diameter but no smaller (i.e., any 
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smaller fragments ejected at higher velocities will escape). 



The second point is that crater measurements by Kirchoff and Schenk (2009), 
which cover a larger region of heavily cratered terrains, show a similar behav- 
ior, namely a decreasing crater spatial density below diameters of a few km, 
and variations in density at larger crater diameters. Some cratered regions are 
at saturation densities, while others are not. It is possible that a combination 
of saturation effects, and a production function that is more steeply sloped at 
larger diameters, leads to the observed crater distributions. 



4-2 The Non- Detection of Craters at the South Pole 

The images provide essentially complete longitudinal coverage of the south 
pole between latitudes of -75° S and -90° S. We did not identify a single crater 
in this region at the limiting image scale of ~ 123 m/pixel, or a completeness 
limit of just over 600 m (i.e., ~ 5 pixels). Based on the lack of observed craters 
in this region, we set an upper limit on crater density. The area of a spherical 
cap is Ac = 2ttR^{1 — sin0), where Rm is the moon's radius and is the 
latitude; thus the area of the region southward of -75° S is 13, 607km^. 

If n is the true crater density (number of craters per square km) and 
is the surface area, then the expected number of craters is < N >= nAc- 
Because primary cratering is a Poisson random process, the probability that 
the observed number of events equals k is: 

P{N = k) = (19) 
where /i is the expected value. In our case, fj, = nAc, and k = 0. Thus the 
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probability becomes: 



P{N = 0) = = e-"^^ (20) 

For a 99% confidence level, or P = 0.01, we derive an upper limit on the true 
crater density of = 3.38 x 10~^ craters/km^. For a diameter bin size from 
Dmin = 600 m to Dinax = 600^2 ~ 849 m, and bin "middle" of D^id = 683 m, 
the corresponding upper limit on the R-value for the south pole of Enceladus 
is: 

Tl X 

Rs, = J. _n = X 10'' (21) 

This value resides within the range of densities measured in the other young 
terrains. 



4-3 The Primary Crater Size-Frequency Distribution 



The primary crater SFD on the young surface of Europa (Bierhaus et al. , 2009 



Schenk et al. , 2004 ) displays an approximately -2 differential slope at crater 



diameters between a few km and a few tens of km, as do younger terrains on 



Ganymede (bright terrain and floors/ejecta blankets of impact basins, Schenk 



et al. (2004)) and Callisto (floors/ejecta blankets of impact basins, Schenk 



et al. (2004)). The behavior of the primary crater SFD at smaller diameters is 



masked by the extensive secondary populations (Bierhaus et al. , 2005). Thus 



Enceladus's young terrains provide an important test of the primary crater 
SFD in the outer solar system at diameters smaller than a few km. 



The crater SFDs on the young terrains of Enceladus have a varying but con- 
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sistent trend of an approximately -2 differential slope. The density varies by 
up to an order of magnitude between regions; the lowest crater density is at 
and near the south pole, while the higher crater densities are at the mid- to 
northern latitudes. The consistency in the crater SFD, within young terrain 
and regardless of crater density, strongly suggests that this shape is in fact 
the production population. 



The south polar plumes deposit material on the surface, and Kirchoff and 



Schenk (2009) find morphological evidence that some craters do have a "soft- 



ened" appearance. However, deposition from the plumes cannot account for 



the shallow small crater SFD we see on Enceladus. Kempf et al. (2010) predict 
maximum deposition rates of about 1 mm/yr, but only for locations within 
100 m of a plume. For more distant locations, beyond 10 km, the deposition 
rate drops off to less than 10^^ mm/yr, and the majority of the area that 
sees plume deposits have deposition rates more like 10~^ mm/yr. A 500 m 
diameter crater with a 100 m depth would take 1 Gyr to disappear by plume 



infill. Meanwhile 500 m diameter craters, using the estimates of Zahnle et al. 



(2003), are expected to form on Enceladus on timescales of 0.01-1 Myr, or at 
least three orders of magnitude faster. Thus we conclude that the small crater 
SFD on Enceladus young terrains (largely) reflects the production population. 



The Enceladus young terrain crater SFD, combined with observations from the 
Galilean satellites, provide preliminary cross-planetary evidence for a ~ —2 
differential-slope, primary-crater population at crater diameters less than a few 
tens of km in the outer Solar System. This argues for an impacting population 
(i.e., comets) that have a similar SFD for projectile diameters in the range from 
< 100 m to several km. 
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4 -4 Mimas 



We measured an image mosaic that includes the large (~ 140 km) Herschel 
crater and surrounding region. Table 5 lists the images, and Figures 22 and 
23 show the mosaic and R-plot, respectively. The measurements are for the 
region outside Herschel crater. Much like the measurements on the heavily 
cratered terrain of Enceladus, the crater SFD for diameters larger than 5 km 
has a roughly -3 differential slope (which is fiat on an R-plot). At smaller 
diameters, the crater density decreases, and the SFD transitions to a roughly 
-2 differential slope. 

Because the measurements are in the region immediately outside the Herschel 
crater, it is not unreasonable to expect that the large crater has affected the 



local small crater population. However, we note that Kirchoff and Schenk 



(2010) measured a large region of Mimas that includes areas not adjacent 
to Herschel, and find a very similar SFD (compare their Figure 2 with our 
Figure 23). 



4-5 Rhea 

We measured the crater SFD in images located near the leading face of Rhea. 
Table 6 lists the images used to make the mosaic seen in Figure 24, and 
Figure 25 is the R-plot of our data. We measured almost 7500 craters in just 
the four images outlined in Figure 24. 

Unlike Mimas and Enceladus, the crater density does not decrease at diameters 
smaller than a few km; indeed. Figure 25 shows an increase in crater density 
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at smaller diameters. This observation matches the prediction of Section 2.3, 
namely that Rhea (along with Tethys and Dione) is sufficiently massive that 
it can retain a measurable secondary crater population. This is explicitly illus- 
trated by an image sequence on Rhea that captures (one of the few examples 
imaged to date on Saturnian satellites) a secondary crater cluster near the 



48-km ray crater Inktomi (Wagner et al. , 2011). 



5 Sesquinaries: Where Art Thou? 

The correlation between predicted and observed secondary crater populations 
on the Saturnian moons is encouraging. What remains puzzling is the absence 
of a measurable signature of sesquinary craters on Enceladus, and especially on 
Mimas. The fragments that make secondary craters follow an inverse mass- 



velocity relationship (e.g., McEwen and Bierhaus (2006)), which partly ex- 



plains their generally steep SFDs. Because sesquinary fragments from moons 
with low escape velocities are presumably the same fragments that make sec- 
ondary craters on moons with higher escape velocities, one would logically 
expect that Mimas, with its multitude of large craters, would be covered in 
a steeply-sloped population of small craters made from sesquinary fragments. 
Yet the crater SFD clearly show a decrease in crater density at smaller sizes 
on an R-plot. Some of the small craters may, in fact, be sesquinaries, and at 
the moment we have no means to identify whether a small crater is primary 
or sesquinary. But the fact remains that the amount of escaped ejecta (see Ta- 
ble 3) per impact is factors of several greater than the impactor mass. Given 



that the sesquinary fragments mostly re-impact their source moon (Alvarellos 



et al. , 2005 ), there should be a high-density small crater population made from 
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sesquinary fragments. What happens to those fragments? 



Zahnle et ah (2008) estimate that Ionian sesquinaries on Europa should dom- 



inate the small crater population on Europa, provided that Melosh's spall 
model can be applied without modification to the 2.5 km/s ejection velocities 



required for spalls to escape lo. Yet observations (Bierhaus et al. , 2005) clearly 



show that most small craters on Europa are Europan secondaries. Zahnle et al. 



(2008) wondered whether the Melosh (e.g., Melosh (1984)) spall model applied 
to the higher velocity fragments that are likely the source of the sesquinary 
population. The spall model predicts a plate-like geometry for the fragments, 
i.e. they are much larger in length and width than they are in height. As 



Zahnle et al. (2008) discussed, it may be unreasonable that such fragments 



survive in that state; more reasonably, they may break up into smaller pieces 
whose mean radius is approximated by the plate thickness, or by some smaller 
size scale intrinsic to the target material (for lo this might be the thickness of 
individual lava flows, presumably of order a meter). 



Schultz and Gault (1985) conducted a number of experiments that explored 



the morphology of craters made by clusters of fragments rather than individual 
fragments. The craters made by fragments that were still relatively compact 
appeared similar to traditional impact craters, albeit shallower (as secondaries 
tend to be). More dispersed fragment clusters made craters whose morphology 
departed from a traditional impact crater. 



An explanation for the missing sesquinary craters may lie in the mechani- 
cal construction of the ejecta fragments, and the greater (much greater, in 
many cases) time between launch and impact for sesquinary craters and sec- 
ondary craters. Adjacent secondary craters are made by fragments with short 
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flight times, less than a few minutes and sometimes less than a minute. Dis- 
tant secondary craters are made by fragments with flight times of up to tens 



of minutes, depending on the moon. However, Alvarellos et al. (2005) show 
that sesquinary fragments in the Saturnian system typically have lifetimes of 
decades - some even survive for 10,000 years. This is between thousands to 
one million orbital periods of the moons, and thus the fragments have orders 
of magnitude more time to disperse than do fragments that make secondary 
craters. Close encounters with the parent moon may provide additional dis- 
persive pulses to weakly bound fragments. The smallest discrete fragment, 
resistant to any further disruption or dispersion, flnally impacts its source 
moon; perhaps its size is sufficiently diminutive that the resulting crater is 
below the imaging resolution of the current Cassini ISS data. For example, 
a 300-m crater on Enceladus, which we would probably see in the highest- 
resolution mosaics (Table 4), would be created by an 85-m diameter fragment 
striking at the moon's escape velocity, assuming a 45° impact angle (Eq. 13). 
For a km-sized crater on Enceladus, which we would see in any of the mo- 
saics, the corresponding impactor diameter is 380 m. These impactor sizes are 
intermediate between the characteristic sizes for rubble and spalls calculated 



by Alvarellos et al. (2005) (see their Figures 2 and 4), so it seems plausible 
that modest fragmentation in orbit could render the sesquinaries unobserv- 



ably small. See Alvarellos et al. (2008) and Zahnle et al. (2008) for further 
discussion of fragment sizes. 



6 Summary 



The following are the key observations and outcomes of this analysis: 
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(1) The mass available to make secondary craters depends upon the relative 
magnitudes of f (the minimum velocity required to make a secondary 
crater) and the target moon's modified escape velocity f^^,, in addition to 
impactor and target mechanical properties (e.g., material strengths and 
density contrasts). Objects for which v^^^ < v^i^ should have no secondary 
crater population. For the current discussion, that regime includes Mimas 
and perhaps Enceladus (depending upon the appropriate value of f mm for 
Enceladus - some variability in the small crater SFD seen on Enceladus' 
young terrains provides preliminary evidence that secondaries may in fact 
form on Enceladus). More broadly, this applies to all small moons and 
most minor planets. 

(2) A low surface gravity also means that primary craters will lack an adja- 
cent secondary crater population, as even low-velocity ejecta will travel 
far from the parent primary, perhaps even significant fractions of the 
body's circumference. 

(3) Rhea, Tethys, and Dione are sufficiently massive that we expect them 
to retain measurable secondary crater populations. Our measurements 



of Rhea (as well as those of Kirchoff and Schenk (2010)) support this 
prediction. However, these moons will lack dense populations of adjacent 
secondaries due to the process described above in (2). 

(4) lapetus should not have significant secondary or sesquinary crater pop- 
ulations. This is a consequence of the low primary impact velocities due 
to lapetus' large distance from Saturn, which overall generate smaller 
primary craters with less ejecta, and also generate slower-moving ejecta. 

(5) The Galilean satellites, and Europa in particular, are in the "sweet spot" 
for secondary crater production. High impact velocities generate large 
craters with lots of fast-moving ejecta, and the moons are massive enough 
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to retain that ejecta to form secondaries. The same is true for the Moon, 
Mars, and Mercury. 
(6) The crater SFD measured on young terrains of Enceladus generally fol- 
lows a -2 differential slope. The similarity between this SFD on Enceladus, 
the primary crater SFD on the young jovian moon Europa, and recent 



observations of young rayed craters on other Saturnian satellites (Schenk 



and Murphy, 2011) provide growing evidence that this is, in fact, the 



production SFD at craters of these sizes, and that comets making these 
craters (comets less than a few hundred meters in diameter) follow a 
similar distribution. A slope of -2 for the nuclei of Jupiter-family comets 



(JFCs) was recently pubhshed by Weiler et al. (2011), but this result 
cannot be taken as a confirmation of our work because (1) there are very 
few determinations of nuclear sizes for sub-km comets; (2) few JFCs have 
been discovered by surveys with well-characterized observational biases; 
and (3) since JFCs have, by definition, been active at some point, they 
have undergone physical evolution, which will tend to flatten their size 
distribution. Triton appears to have a signiflcantly steeper crater SFD 



(Stern and McKinnon, 2000 McKinnon and Singer, 2010), but its huge 
apex-antapex asymmetry seems inconsistent with heliocentric impactors 



(Schenk and Zahnle, 2007), thus making the link between crater sizes and 
the SFD of the impactors a thorny one. 
(7) The fate of sesquinary ejecta remains elusive. There is little doubt that 
large impacts on small satellites send signiflcant amounts of mass into or- 
bit, mass that would otherwise make secondary craters on larger moons. 
However, unlike secondaries - which leave a distinct imprint on small 
crater SFDs by adding steep slopes - sesquinaries do not express them- 
selves so ostentatiously. Heavily cratered regions on Enceladus, and espe- 
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cially Mimas, should have small crater populations near or at saturation 
density due to the effect of sesquinaries, and yet both moons display a 
distinct decrease in small crater density. One explanation may be that the 
fragments responsible for secondary craters are not yet sub-divided into 
their smallest discrete components. The significantly longer time spent as 
sesquinary fragments allows the ejecta fragments to further separate into 
their smallest discrete components, which form craters too small to be re- 
solved by current image data. On the other hand, Cassini discovered three 
small moons - Pallene, Methone, and Anthe - with diameters of ~ 5, 3, 



and 2 km, respectively, orbiting between Mimas and Enceladus. |Alvarel- 



los et al. (2005) speculated that these moons might be spalls launched in 



an event like the Herschel impact. Thus there is some reason to believe 
that large fragments can indeed be launched (or accrete in orbit around 
Saturn) and survive for long periods of time. 
(8) The bulk of the crater SFDs across the Saturnian satellites may be ex- 
plained by a single impacting population. The variation in impact veloc- 
ity and surface gravities across the moons means that a single impacting 
population will generate different primary crater SFDs on each satellite. 
For a given sized impactor, the variation in primary crater size is followed 
by variation in ejecta mass and ejecta speeds available to make secondary 
and (probably) sesquinary craters. For example, on lapetus a 1 km comet 
makes an approximately 8 km transient crater with ~ 2 x 10^^ kg available 
to make secondaries, while the same impactor makes a 17 km transient 
crater on Mimas, with no mass available to make secondaries. The mass 
available to make secondaries will travel different distances across the 
moons. The superposition of the varying primary crater, and resulting 
secondary (and sesquinary) crater distributions may explain the crater 



37 



SFDs seen on the satellites. 
(9) Finally, we propose an update to the Voyager-era interpretation of the 
Saturnian cratering record: Population I is likely dominated by helio- 
centric comets, and appears on all Saturnian satellites as the source of 
craters above several km diameter; Population 11 is a result of secondary 
(and perhaps sesquinary) craters, with significantly varying signatures 
between the satellites, due to differences in primary impact velocities, 
surface gravities, and escape speeds. 
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Moon vj [km/s] fesc [m/s] v^^^ [m/s] g [m./s^] 

Mimas 27 159 130 0.064 

Enceladus 24 239 209 0.11 

Tethys 21 393 345 0.145 

Dione 19 510 467 0.231 

Rhea 16 635 559 0.285 

lapetus 6.1 572 566 0.223 

Europa 26 2000 1888 1.31 

Ganymede 20 2741 2629 1.43 

Callisto 15 2440 2382 1.235 

Table 1 

Table of satellite properties used in our calculations. ^ Impact velocities from Zahnle 

et al. (2003). is the modified escape velocity, see Equation 4. 



47 



Moon D [km] Mtot Mbik /bik M^ec /sec Mi.5 /: 



1.5 



Mimas 17.3 17.485 x 10^^ 17.294 x lO^^ o.989 0.000 0.000 19.065 x 10^^ 0.011 

Enceladus 14.9 11.267 x 10^^ 11.134 x lO^^ 0.988 4.850 x 10^^ 0.004 8.430 x lO" 0.007 

Tethys 13.5 8.307 x lO^^ §.196 x 10^^ o.987 7.547 x 10^^ 0.009 3.498 x lO" 0.004 

Dione 11.9 5.717 x lO^^ 5 g2l x lO^^ 0.983 7.616 x lO" 0.013 2.005 x lO" 0.004 

Rhea 10.7 4.191 x lO^^ 4 115 ^ lo^s 0.982 6.455 x lO^^ 0.015 1.135 x 10^^ 0.003 

lapetus 7.8 1.630 x 10^=^ 1.610 x lO^^ 0.988 1.685 x 10^^ 0.010 0.321 x lO^^ 0.002 

Europa 9.7 3.071 x 10^^ 2.922 x 10^^ q 952 14.351 x 10^^ 0.047 0.450 x 10^^ 0.001 

Ganymede 8.6 2.178 x 10^^ 2.075 x 10^^ q 953 iq 123 x 10^^ 0.047 0.198 x 10^^ 0.001 

Calhsto 8.0 1.712 x 10^^ 1.643 x 10^^ q 959 6.790 x 10^^ 0.040 0.153 x 10^^ 0.001 
Table 2 

The total and fractional masses for the primary crater and resulting ejecta when 
■Wmin = 150 m/s. The diameter listed in the second column is the transient crater 
diameter predicted by Equation 13, and not the final crater diameter. The subscript 
hlk refers to the ejecta blanket, the subscript sec refers to secondary craters, and the 
subscript 1.5 refers to sesquinaries. This case of Umin minimizes Mbik, and yet almost 
all the mass is still moving slower than Vmin, i-e., most the mass is still in Mbik- All 
masses are in kg. The fractional masses are /bik = Mbik/Mtot, /sec = -^sec/-^tot) 
and /1.5 = Mi.s/Mtof 
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^^min = 150 m/s 




^min = 250 m/s 
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Tethys 


2.4 


1.1 
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1.1 


Dione 


2.4 
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0.6 








u.o 


u.^ 


lapetus 


0.5 


0.1 


0.2 
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Europa 


4.6 


0.1 


2.2 


0.1 


Ganymede 


3.2 


0.1 


1.6 


0.1 


Callisto 


2.2 


< 0.1 


1.0 


< 0.1 



Table 3 

Msec SLnd M1.5 in terms of the impactor mass, rrii, when v„im = 150 m/s and 
Vmin = 250 m/s. The horizontal lines in the table divide the moons into groups 
based on the amount of ejecta available for secondaries and sesquinaries. See text 
for discussion. 
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Enceladus Mosaic Image Scale [m/pix] 

ISS_003BN_LIMTOP004_PRIMB N1487300107_l 148 

N148730028B_1 141 

N1487300482_l 133 

N1487300648-1 126 

N1487300854_l 118 



ISS_004BN_RBGEO002_PRIME N1489047900_2 166 

N1489048083_2 158 

N14890482BB_2 151 

N1489048757_2 131 

N1489048931_2 124 

N148904910B_2 117 



ISS_004BN_RBGEO002_PRIME N1489049969-2 82 

N1489050144_2 76 

N1489050320_2 68 

N148905047B_2 62 

N14890506B1_2 56 



ISS_011EN_MORPH002_PRIME N1500062131_l 80 

N1500062262_l 73 
N1500062382_l 67 



ISS_011EN_N9COL001_PRIME N1500061132_l 129 

N1500061253-1 123 

N1500061390_l 116 

N1500061612_l 110 

N1500061634_l 104 

N1500061771_l 97 

N1500061892_l 91 

Table 4 

The mosaics and individual images we used to make crater SFD measurements on 
Enceladus. 
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Mimas Mosaic Image Scale [m/pix] 

ISS_126MI_GEOLOG001_PRIME N1644777693_l 93 

N1644777828_l 97 

N1644777993_l 102 

N1644778141_l 106 

N1644778308_l 110 

N1644778455_l 115 

N1644778567_l 118 

Table 5 

The images used to measure the crater size-frequency distribution on Mimas. 
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Rhea Mosaic 



Image Scale [m/pix] 



ISS_121RH_REGMAP001_PRIME 



N1637518901_l 147 

N1637519058_l 148 

N1637519176_l 148 

N1637519283_l 149 

N1637519392_l 151 

N1637519501_l 152 

N1637519610_l 154 

N1637519768_l 156 

N1637519875_l 159 



Table 6 

The list of images used to measure a representative population of Rhea's crater 
size-frequency distribution. 
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1. Primary 2. Spmndary 3. Sp^quinary 




Fig. 1. A simple schematic of the three types of impact craters considered in this 
paper. (1) Primary craters are caused by the direct impact of a comet or an aster- 
oid. (2) Secondary craters are caused by ejecta from a primary crater that travels 
ballistically to some distance away, and which impacts the surface with sufficient 
speed to form a crater. (3) Sesquinary craters are caused by ejecta that initially 
escapes the satellite, and go into temporary orbit around the planet. Some time 
later the material the material impacts a satellite (usually the moon from which the 
material was ejected) to form a crater. 
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Fig. 2. A schematic of the amount of mass ejected faster than a given velocity. There 
is an inverse relationship between mass and velocity, i.e. less mass is ejected at higher 
velocities. "Umin is the minimum velocity required to make a secondary, while t^esc 
is the escape velocity for a body. Mass ejected faster than the escape velocity is 
the same as the mass available to make sesquinary craters, or M{v > v^^) = Mi. 5. 

is the modified escape velocity, see text for discussion.) Mass ejected faster 
than Vmin but slower than fesc is the mass available to make secondary craters. Msec 
(unless Vesc < Vmin, in which case no mass is available to make secondaries). Mass 
ejected slower than fmin is available to make an ejecta blanket, Mbik- 
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Fig. 3. A Galileo image of Rhiannon, which is a ~ 15 km diameter primary crater 
on Europa. The nearest distinct secondary craters are about 17 km from a point 

half-way between Rhiannon's center and the crater rim. 
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Fig. 4. A Galileo image of Tyre, which is a 45 to 50 km diameter primary crater on 
Europa, depending on which feature is the true crater rim. The secondary crater 
population begins about 50 km from the central portion of the crater. The black 
bar is a data gap. 
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Fig. 5. A plot of Msec (top panel) and /sec (bottom panel), the absolute and fractional 
ejected mass available to make secondary craters for a 1 km cometary impactor, for 
"ymin = 150 m/s. The Saturnian satellites are open circles, and the Galilean satellites 
are closed circles. The mass available for secondaries on Mimas is zero for this value 
of Vrain, i-G. none of the ejecta retained by the satellite are moving fast enough 
to make secondary craters. The Galilean satellites, and Europa in particular, have 
much more mass available for secondary craters. 
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Fig. 6. A plot of Msec (top panel) and /sec (bottom panel), the absolute and fractional 
ejected mass available to make secondary craters for a 1 km cometary impactor, for 
Vmin = 250 m/s. The Saturnian satellites are open circles, and the Galilean satellites 
are closed circles. The mass available for secondaries on both Mimas and Enceladus 
is zero for this value of Vmm, i-e- none of the ejecta retained by the satellites are 
moving fast enough to make secondary craters. 
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Fig. 7. A plot of Mi.5 (top panel) and /1.5 (bottom panel), the absolute and fractional 
ejected mass available to make sesquinary craters, for a 1 km cometary impactor. 
The Saturnian satellites are represented by open circles, while the Galilean satellites 
are shown as closed circles. M1.5 is independent of Vmin- Mimas and Enceladus have 
much more mass available to make sesquinary craters than the other Saturnian 
satellites or the Galilean satellites. 
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Fig. 8. The top plot is the non-planar ballistic range in km on each of the moons for 
ejecta launched at 150 m/s and a 45° angle, as a function of the surface gravities 
of the satellites. The bottom plot is that range, normalized to the circumference of 
each satellite. Mimas does not appear because ejecta launched at 150 m/s escapes 
that moon. The horizontal axes are the same in the two plots. See text for further 
discussion. 
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Fig. 9. The same as Figure 8, but for ejecta launched at 250 m/s. The bottom plot 
is that range, normalized to the circumference of each satellite. Neither Mimas nor 
Enceladus appear because ejecta launched at 250 m/s escapes these moons. The 
horizontal axes are the same in the two plots. See text for further discussion. 
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Fig. 10. (A) is the full frame, unprocessed Cassini image N1500062262_l, which is one 
of the images of the ISS_011EN_MORPH002_PRIME mosaic. (See Figure 18.) The 
image is complete on the left-hand side, but every-other line on the right-hand side 
ends prematurely due to the data compression algorithm. The termination point of 
every-other line is different for every line pair, because different scene content within 
each line pair will compress with different efficiencies. The white box outlined in 
the upper-right is a 250 x 250 pixel portion of the image. (B) is simply a slightly 
magnified portion of the region within the 250 x 250 white box, better illustrating 
the every-other-line data gaps. (C) is the gap-filled version of (B), accomplished 
by using the ISIS routine lowpass. Essentially, each blank pixel is populated by an 
intensity value that is the average of the intensity values in the pixels immediately 
above and below the blank pixel. 
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Fig. 11. The outlines of the Enceladus mosaics we measured. The latitude lines are 
in increments of 15°, and the longitude lines are in increments of 30° and are in 
positive-east longitude. The measured regions span latitudes from the south pole to 
a little north of the equator, and span a variety of longitudes. The shadowing on 
the sphere is arbitrary, and does not reflect the terminator in any of the measured 
image sequences. 
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Fig. 12. Some of the higher-resolution images of the Enceladus 
ISS_004EN_REGEO002_PRIME mosaic. The average mosaic scale is 69 m/pix, 
and the maximum scale is 83 m/pix. We divided the crater measurements into the 
young, fractured terrain (outlined in blue) and the heavily cratered terrain. 
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Fig. 13. The R-plot for the regions seen in Figure 12. The measurements within 
the young, fractured terrain are the open blue symbols, while the measurements 
within the older, cratered terrain are the filled symbols. We measured 15 craters in 
the fractured terrain, and 2375 craters in the older terrain. The fractured terrain 
displays a clear ~ — 2 differential slope. The cratered terrain has a ~ —3 differential 
slope above a few km, which transitions to a ~ — 2 slope at smaller diameters. See 
text for discussion. 
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Fig. 14. Some of the lower-resolution images of the Enceladus 
ISS_004EN_REGEO002_PRIME mosaic. The average mosaic scale is 132 m/pix, 
and the maximum scale is 173 m/pix. We divided the crater measurements into 
two regions of young, fractured terrain (outlined in purple and green) and the 
heavily cratered terrain. 
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Fig. 15. The R-plot for the regions seen in Figure 14. The measurements within the 
young, fractured terrain are the purple and green symbols, while the measurements 
within the older, cratered terrain are the black symbols. We measured 153 craters in 
the fractured region outhned in purple, 26 craters in the fractured region outlined 
in green, and 5181 craters in the cratered region. Though the R values show a less 
clear linear trend than the young terrain measurements seen in Figure 13, the young 
terrain SFDs here too follow a ~ —2 differential slope. The cratered terrain has a 
~ —3 differential slope above a few km, which transitions to a ~ —2 slope at smaller 
diameters. See text for discussion. 
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Fig. 16. Some of the images of the Enceladus ISS_003EN_LIMTOP004_PRIME mo- 
saic. The maximum image scale is 150 m/pix. We divided the crater measurements 
into two regions of young, fractured terrain (outhned in blue) and the heavily 
cratered terrain. 
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Fig. 17. The R-plot for the regions seen in Figure 16. The measurements within 
the young, fractured terrain are the blue symbols, while the measurements within 
the older, cratered terrain are the black symbols. We measured 48 craters in the 
fractured terrain, and 1840 craters in the cratered terrain. Again, the young ter- 
rain SFDs here follow a ~ —2 differential slope. The cratered terrain has a ~ —3 
differential slope above a few km. See text for discussion. 
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Fig. 18. The Enceladus ISS_011EN_MORPH002_PRIME mosaic. The average mo- 
saic scale is 67 m/pix, and the maximum image scale is 80 m/pix. We divided the 
crater measurements into two regions of young, fractured terrain (outlined in blue 
and green) and the heavily cratered terrain (outlined in purple). 
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Fig. 19. The R-plot for the regions seen in Figure 18. The measurements within the 
young, fractured terrain are the blue and green symbols, while the measurements 
within the older, cratered terrain are the purple symbols. We measured 4 craters 
in the fractured region outlined in green, 15 craters measured in the fractured area 
outlined in blue, and 53 craters in the cratered terrain outlined in purple. The young 
terrain SFDs here follow a ~ —2 differential slope. The cratered terrain area is too 
small to contain large craters with a spread of diameters, and so does not have a 
distinct slope for diameters above several km. However, like the heavily cratered 
regions seen in previous figures, the crater density decreases at smaller diameters. 
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Fig. 20. The Enceladus ISS_011EN_N9COL001_PRIME mosaic. The average mo- 
saic scale is 110 m/pix, and the maximum image scale is 129 m/pix. We show a 
latitude/longitude grid on this mosaic to emphasize the extent of the young terrain 
imaged in this mosaic, and the excellent coverage of the south pole itself. The 180 
label refers to the longitude of that line. For this mosaic, we measured only the 
young terrain. 
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Fig. 21. The R-plot for the young terrain seen in Figure 20. We measured 22 craters 
in this region. The ^ —2 differential slope is similar with other measurements within 
the fractured terrain, but at a lower density (lower R- value), consistent with the 
expectation that the area around the south pole is the youngest region on Enceladus. 
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Fig. 22. The Cassini ISS_126MLGEOLOG001_PRIME mosaic of Mimas, with an 
average image scale of about 105 m/pix. The large crater in this mosaic is Herschel 
crater, about 135 km diameter. The measurements we report are for the area outside 
Herschel crater (in between the blue and black outlines). 
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Fig. 23. An R-plot for our measurements of the area shown in Figure 22; we mea- 
sured 6581 craters in that region. The density decrease at crater diameters less 
than about 10 km is not a completeness effect, because the average mosaic scale 
is about 105 m/pix. The decrease in small crater density supports our calculation 
that secondary craters should be rare to non-existent on Mimas. 
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Fig. 24. A subset of the Cassini ISS_121RH_REGMAP001_PRIME Rhea mosaic, av- 
erage mosaic scale of about 159 m/pix. The presence of a higher-density small crater 
population supports our prediction that Rhea is sufficiently massive to support the 
production of a secondary crater population. 
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Fig. 25. An R-plot for our measurements of the area shown in Figure 24; we measured 
7495 craters in that region. Unhke Enceladus or Mimas, there is not a crater density 
decrease at smaller diameters. 



77 



